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We study theoretically dynamics in a Josephson junction coupled to a mechanical resonator 
looking at the signatures of the resonance in d.c. electrical response of the junction. Such a system 
can be realized experimentally as a suspended ultra-clean carbon nanotube brought in contact with 
two superconducting leads. A nearby gate electrode can be used to tune the junction parameters 
and to excite mechanical motion. We augment theoretical estimations with the values of setup 
parameters measured in the samples fabricated. 

We show that charging effects in the junction give rise to a mechanical force that depends on the 
superconducting phase difference. The force can excite the resonant mode provided the supercon- 
ducting current in the junction has oscillating components with a frequency matching the resonant 
frequency of the mechanical resonator. We develop a model that encompasses the coupling of elec- 
trical and mechanical dynamics. We compute the mechanical response (the effect of mechanical 
motion) in the regime of phase bias and d.c. voltage bias. We thoroughly investigate the regime of 
combined a.c. and d.c. bias where Shapiro steps are developed and reveal several distinct regimes 
characteristic for this effect. Our results can be immediately applied in the context of experimental 
detection of the mechanical motion in realistic superconducting nano-mechanical devices. 

PACS numbers: 85.85. +j, 74.45. +c, 73.23.Hk 



I. INTRODUCTION 

Nanoscale electromechanical systems (NEMS) convert 
small amplitude mechanical motion into measurable elec- 
trical currents [T]. Devices based on NEMS have found 
applications as sensitive detectors of mass force [3] 
and electrical charge j4j. Considerable research efforts 
have been dedicated to improving detection sensitivity 
by fabricating devices with higher resonance frequencies, 
lower damping rates (high quality factors) and larger cou- 
pling between electrical and mechanical degrees of free- 
dom. 

The problem of detecting the quantum state of a 
macroscopic mechanical resonator gave rise to several 
measuring schemes, proposed [5J [B] as well as realized 
[7]. Continuous improvements in device fabrication in 
combination with new techniques for cooling mechanical 
motion [51 [S] have pushed the sensitivity threshold to the 
quantum limit [7]. The use of superconducting devices, 
in particular, superconducting qubits to detect and con- 
trol the mechanical motion is in focus of modern research 
[3 [TU] . It gives rise to a growing interest in techniques 
of coupling NEMS to superconducting circuits. 

Superconducting nano-devices frequently use Coulomb 
blockade that makes their transport properties sensitive 
to the gate voltages. The same gate voltage can be used 
to excite the mechanical motion which is detected from 
the change of d.c. transport properties of the device [T31- 
[T3] . Without superconductors, this scheme has been suc- 
cessfully realized for a metallic single-electron transistor 
[17j and for a Coulomb-blockaded quantum dot in an 
ultra-clean carbon nanotube (CNT) [T5] . The results re- 
vealed high resonance frequencies, reaching gigahertzs, 
and unprecedented quality factors of the order of 10^. 
These devices can be made superconducting by connect- 



ing them to superconducting leads and providing suffi- 
ciently large coupling between the states of the lead and 
device. We have successfully realized Josephson junc- 
tions based on the ultra-clean CNT. The supercurrent 
observed demonstrates a pronounced gate-voltage sensi- 
tivity that indicates a well-developed Coulomb blockade 

A very interesting proposal that combines Josephson 
dynamics and mechanical resonator has been recently put 
forward by Gothenburg collaboration [20] • The authors 
consider an ideal ballistic CNT between two supercon- 
ducting leads biased at voltage V. Owing to Josephson 
relation, the current in the nanotube oscillates at fre- 
quency Ljj = 2eV/h. The authors notice that in external 
magnetic field this gives rise to an oscillating Lorentz 
force. If the frequency matches the frequency of the me- 
chanical resonator, the force excites mechanical motion 
which rectifies the Josephson current enabling the ob- 
servation of the effect in d.c. electric response of the 
junction. One would observe a narrow current peak in 
I —V characteristics of the device. The same mechanism 
is responsible for Fiske steps [23] '■ the difference is that in 
Fiske experiments the resonance is electrical rather than 
mechanical. 

This provides us motivation for the present theoretical 
study where we address superconducting NEMS where a 
mechanical resonator is integrated with a superconduct- 
ing circuit element, a Josephson junction. The focus of 
this work is to describe the techniques of driving and de- 
tecting mechanical motion using superconducting current 
and to explore the coupled dynamics of the oscillator dis- 
placement and superconducting phase difference. While 
the study is general for any type of oscillating nanowire, 
for the purpose of illustration and for concrete estima- 
tions we refer to one of the more successful NEMS de- 
vices: a suspended metallic CNT resonator connected to 
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the metallic leads. We have fabricated and studied such 
devices. 

Studies of CNT Josephson junctions have shown that 
their Josephson energy can be modulated by the gate- 
induced charge on the CNT. In this paper, we have shown 
that this gives rise to a different and generally more im- 
portant mechanism of mechanical driving than that con- 
sidered in ^20]. We consider all possible non-linearities 
in coupled Josephson- mechanical dynamics. We have 
shown that in our situation the most important one is 
the mechanical non-linearity. We provide detailed esti- 
mations of the displacement, force and electrical current 
scales involved. It is our conclusion that the mechanical 
response in our devices should be small modifying the 
current on the scale of 10~^ of the critical current (In this 
case, the critical current is estimated from the product 
of the junction conductance and superconducting energy 
gap and can exceed the experimentally measured switch- 
ing current by two orders of magnitude) . 

In the present paper, we have investigated in detail the 
mechanical response under conditions of phase bias and 
d.c. voltage bias. We dedicated special attention to the 
dynamics in the presence of external a.c. drive, in the 
regime where Josephson junction exhibits well-developed 
Shapiro steps [21] . One of the motivations of the research 
on Shapiro steps was the better synchronization condi- 
tions in comparison with d.c. voltage bias. The reason for 
this is that the big quality factor Q of the nanomechani- 
cal resonance results in a narrow Fiske-like current peak. 
Its width in voltage can be estimated as 5V ~ V/Q. 
This imposes a severe limitation on voltage noise Sy'- 
to resolve the peak one must achieve Sv ^ {e/h)V/Q. 
This may be challenging under realistic experimental cir- 
cumstances. There is a way out: the voltage can be 
synchronized with the frequency of external irradiation. 
This effectively reduces the noise. From the other hand, 
the external irradiation can excite the mechanical motion 
by itself. The resulting complicated dynamics should be 
augmented with non-linear effects. 

In this paper we present our theoretical results con- 
cerning the mechanical response manifested in d.c. elec- 
tric response of the junction, in particular, the modifi- 
cation of the width and position of Shapiro steps, and 
analyze a variety of distinct regimes that differ in rela- 
tive and absolute magnitudes of driving forces and reso- 
nant conditions, and are manifested in distinct lineshapes 
of the resonant response. Our preliminary experimental 
results show corresponding features. They will be pre- 
sented elsewhere |25j upon completion of detailed analy- 
sis and comparison with our theoretical findings. 

The paper is organized as follows. 

In Section |ll] we describe the setup concentrating sep- 
arately on the electrical and mechanical properties in 
the corresponding Subsections and introduce the nota- 
tions. In Section III we analyze in detail the mecha- 
nisms of coupling between mechanical and electrical de- 
grees of freedom. We consider electrostatic energy, derive 
and explicate the concept of Josephson mechanical force. 



We present a rather involved analysis of competing non- 
linearities and conclude that for our devices the mechan- 
ical non-linearities dominate and the Josephson force can 
be strong enough to induce the non-linear mechanical re- 
sponse. We give the workflow we use to compute the 
mechanical response. In the end of the Section, we spec- 
ify a set of concrete parameters based on experimentally 
measured values. In Section|TV]we address the phase bias 
conditions. We reveal that the phase- dependent shift of 
the resonant frequency can be quite noticeable and dis- 
cuss Lorcntz-like and Fano-like frequency dependences 
of the mechanical response. In Section |V] we discuss the 
Fiske-type mechanical response at Josephson frequency 
matching the resonance frequency of the mechanical res- 
onator, or an integer fraction of this frequency by higher 
harmonics, and give the estimations of the effect. We 
shortly discuss the parametric excitation. In Section VI 
we study the mechanical response at the Shapiro steps 
in the regime where the a.c. driving frequency matches 
the resonant frequency and present the mechanisms and 



peculiarities of this response. In Section VII we consider 
the non-resonant driving that appears to efficiently ex- 
cite the mechanical motion in the regime of Shapiro steps. 



Our concluding remarks are presented in Section VIII 



II. THE SETUP 

The setup under consideration is sketched in Figure [l] 
where we concentrate on a case where both Josephson 
junction and mechanical resonator are realized with the 
same single CNT. Even in this case, the coupling between 
mechanical and electrical degrees of freedom is relatively 
weak. This permits us to describe the electrical and me- 
chanical aspects of the setup separately. We provide such 
description in this Section, while in the next Section we 
concentrate on the coupling. 



A. Electrical Setup 

We consider a conducting link between two supercon- 
ducting leads (the CNT in Fig. 1). In general, the cur- 
rent fiowing in this junction is a complicated non-linear 
and time-delayed function of superconducting phase dif- 
ference between the leads. However we assume that in 
the relevant frequency range the current response of the 
junction is superconducting and instant. 

The junction is included into an external electric cir- 
cuit and the voltage drop at the junction is related to the 
time derivative of Josephson phase, — 2eV/h. In gen- 
eral, a circuit that connects the leads can be described 
by a complex frequency-dependent impedance Zi,(uj) in 
series with a voltage source Vf,- We typically assume that 
Ze by far exceeds the typical junction impedance at low 
frequencies while at frequency scale of Josephson genera- 
tion frequency lu ~ eV/h, Zf^{uj) is negligible in compari- 
son with the junction response. In this case, the junction 
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is current-biased at low frequencies with /(, = Vb/Ze{0) 
and voltage-biased at Josephson frequencies. While this 
scheme looks different from the traditional RSJ model 
where the external impedance is connected in parallel 
and the junction is current-biased, it is equivalent to a 
generalized RSJ upon transforming the impedance and 
the voltage source. For instance, the linear part of pos- 
sible quasiparticle response of the junction can be incor- 
porated into Ze{i-u). 

In addition, the junction is affected by the gate elec- 
trode biased by voltage source Vg. The bias and gate 
circuits are disconnected at zero frequency. At finite fre- 
quency, there is a cross-talk between the circuits which is 
difficult to eliminate or even characterize in realistic ex- 
perimental circumstances. We account for that by corre- 
lating a. c. parts of the voltage sources Vf,,g- For instance, 
if the gate voltage consists of a d.c. part and a harmonic 
signal at frequency ft, Vg{t) = Vgo + Vg cos (fit -I- x), the 
bias voltage source should also oscillate at the same fre- 
quency, Vb{t) — Vbo + VbCos{nt). The ratio of two a.c. 
amplitudes and their mutual phase shift x is determined 
by details of the crosstalk. We will show below that the 
interference of these two a.c. signals may strongly affect 
the d.c. currents in the junction. 

The superconducting current is determined by the in- 
stant phase difference, I{t) = I{(p{t)). In this case, it is 
related to the Josephson energy Ej of the junction, 

Is = {2e/h)dEj{ip)/d^ . (I) 

It is essential for us that the Josephson energy is not 
only a function of phase difference but also depends on 
the gate voltage through the charge q = CgVg induced in 
the resonator, Ej — Ej{q,ip). 

For a nanotube device, the origin of this charge sen- 
sitivity is (weak) Coulomb blockade of electrons in the 
middle of the nanotube. The nanotube can be viewed 
as two junctions in series, those being formed at contact 
with metallic leads. If the conductance of the junctions 
is smaller or comparable with the conductance quantum 
Gq = e^/(7r/i) « 7.75 x I0~^ fi"^ Coulomb interactions 
become important and set a quasiperiodic dependence of 
Josephson energy on q with a period 2e. This corresponds 
to charge quantization in the middle of the nanotube. We 
routinely observe the quasi-periodic modulation of super- 
conducting currents in fabricated nanotube devices. The 
modulation can be tuned by changing the gate voltage 
at scale of g ~ 10 — lOOe from values of the order of 1 
to several per cent. Big modulation and well-developed 
Coulomb blockade require big junction resistances, this 
strongly suppresses the superconducting current. It is 
therefore advantageous to have intermediate resistances 
R ~ Gq^- At i? = 5 fcO we typically observe 30% mod- 
ulation. 

The superconducting current is a periodic function of 
the phase Is{^) = Is{'-P + 27r) and therefore can be ex- 
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FIG. 1: (Color online.) The setup. The sketch presents the 
mechanical resonator realized as a CNT suspended over two 
superconducting leads isolated from the back gate electrode. 
The CNT center is displaced in y-direction by an electrostatic 
force produced by the gate voltage. The superconducting 
leads are parts of an electrical circuit characterized by an 
impedance Z^- The setup can be biased by either voltage or 
current source. 

pandcd in harmonics as |21j : 

oo 

Is{q,(p) ^ h (q) sm{(p) + ^ /„ {q) sm{n(p) , (2) 

If one neglects all harmonics except the first one, Ji 
gives the critical current of the Josephson junction. We 
will typically assume this, and will mention the effect of 
higher harmonics only if it is crucial. 

B. Mechanical Setup 

Mechanical resonators can be realized in a variety of 
ways |11) . In many cases the adequate description of 
the resonator can be achieved with a minimum model 
that accounts for excitations of a single resonator mode, 
neglecting coupling to any other modes. The minimum 
model is given by the following equation of motion for a 
displacement variable y: 

y + Ty + u:ly~ay^~(3v^^F[t)/M . (3) 

Here Fit) is the time-dependent driving force, M is the 
effective mass corresponding to the mode, cjo stands for 
the resonant frequency, F ^ is the damping rate, 
and P is the parameter describing the leading cubic non- 
linearity |12j . The cubic non-linearity provides the im- 
portant restriction on the magnitude of the displacement 
at resonant frequency as a reaction on resonant force. We 
also keep the second-order non-linearity a. Although it 
is not important in analysis of the reaction at resonant 
force, it describes the shift of the resonant frequency due 
to constant force. 

Our preferable realization of mechanical resonator is 
a suspended ultra-clean CNT [TBHTS] that demonstrates 
best quality factors observed so far(Q = wq/F ~ 10^). In 
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this Subsection we review the parameters of the minimal 
model for this realization. In the setup shown in Figure 
[T] the nanotube displacement from equilibrium position 
and the driving force are in the y-direction towards the 
gate, that is, perpendicular to the nanotube axis. The 
mechanical variable y{t) is the displacement of the mid- 
point of the nanotube. 

In the case of a CNT, the adequate model of me- 
chanical properties involves a suspended thin cylindrical 
rod clamped at both ends where the nanotube touches 
the metal leads. The parameters are the rod length 
L, the cylinder radius r, and the tube cross section 
area S. In our experiments, L ~ 0.3..0.5 fim, r ~ 
1 nm, and S — 2iTra ~ 2.1 nrri^ for a single-wall nan- 
otube, a ~ 0.34 nm being the layer spacing in graphite. 
The relevant elastomechanical material constants, car- 
bon Young's modulus E and graphite density p are esti- 
mated as ~ 10^2 [16] and p ~ 2.2 g/cm^. The 
bending modes of the rod and their complete dynamics 
are described by the Euler-BernouUi equation of motion 

nam]. 

We concentrate on the lowest frequency bending mode 
that has no nodes in the rod and therefore is easy to 
excite. The resonant frequency can be tuned by "tight- 
ening" the tube, that is, changing the elastic tension. 
This is achieved by applying a sufhciently big d.c. gate 
voltage Vgo- The resulting electrostatic force strives to 
elongate the nanotube, thus producing the tension. In 
such a way, the resonant frequency can be increased by 
a factor of three in comparison with that of the "loose" 
nanotube. For estimations, we concentrate on the case 
of untightened rod. In this case, the resonance fre- 
quency corresponding to the lowest CNT bending mode 
can be estimated in terms of the bending spring constant 
and the carbon mass density ujq ~ 22 Ay/ EI / pSL^'^ 
|18j . where / — Sr^ is the principal moment of iner- 
tia. In our devices of length L = 0.3..0.5 /j,to, the fre- 
quency is wo/27r ~ 0.30. .0.84 GHz, similar to frequen- 
cies reported in [TS]. The effective force is evaluated 
using the eigenfunction of the mode ^{x) = y{x,t)/y{t), 

F = dxf{x)^{x), f{x) being the force per unit length. 
For electrostatic forces, an ad-hoc assumption is that the 
force distribution is uniform, so the total force is Ft = fL. 
In this case, F ~ 0.53i^t. The effective mass is given by 
M = pS J^dx^^ix), M ~ OAlpSL ~ 4.1. .6.8 x 10-^2 
kg. The cubic non-linearity originates from the tension 
produced by the nanotube displacement, the correspond- 
ing parameter can be estimated as /3 ~ 40 ES/ML^ ~ 
Wp/r^ ~ 1.8. .5.5 GHz^nm"^, assuming uniform distribu- 
tion of force along the length of the rod. The second- 
order non-linearity a vanishes for untightened straight 
rod for symmetry reasons. It however appears if the rod 
is tightened such that its frequency change with respect 
to the untightened value ujq is of the order of ujq. In this 
case, the non-linearity is obtained as a = 3/3j/o, Vo being 
the equilibrium displacement induced by the tightening. 

If F{t) oscillates at frequency uj close to the resonant 
frequency, Eq. ^ can be solved in resonant approxima- 
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FIG. 2: (Color online.) Left panel: Real part of the complex 
displacement amplitude Re y versus detuning. Right panel; 
Imaginary part of the complex displacement amplitude Im y 
versus detuning. The three curves in each panel correspond 
to F/F, = 0.5, 1,2. 



tion for the complex amplitude y: 
F -1 



y = 



2MLJoiy + ir/2 + {p'/2ujo)\y^ 



(4) 



with F being the complex force amplitude. Here we 
introduce the detuning v = u — ujq implying that 
W\ ^ <^o- We also introduce the Duffing parameter 
/?' = P + a^/ujQ ~ oj^/r^ corresponding to the amplitude- 
dependent frequency shift. In our experiments, we esti- 
mate /3' ~ 3.6.. 11 GHz^nm-2. 

We will re-write Eq. Q in dimensionless form intro- 
ducing a critical amplitude yc, yc = \/ I P' ■ At this 
amplitude scale, the response of the resonator becomes a 
two-valued function of detuning (see Fig. |2| For a CNT, 



it can be estimated as y^ 



- IQ which corresponds for 



our experiments to yc — 3.2 •pm. The driving force corre- 
sponding to 2/c is Fc = Mfi'y\ = Moj^ydQ- We estimate 
it for a CNT F^ ~ 10^ ES{r/LfQ~^/'^ ~ 1.2 x IQ-^^ N. 
The dimensionless form of Eq. ^ is: 



y, ir'F, 



R{a,b) = 



2a+\b\^\R{a,b)\^ +i 



(5) 



Here we have introduced a dimensionless complex re- 
sponse function R(y,F). In the linear regime \F\ <C F^ 
its dependence on F can be neglected: i?(a, b) — (2a -I- 

Figure [2] shows the real and imaginary parts of y as a 
function of detuning for three values of the driving force 
amplitude that correspond to quasi-linear, critical and 
bistable regimes. 



III. COUPLING AND NON-LINEARITIES 

In this Section, we analyze the coupling between me- 
chanical and electrical degrees of freedom, with empha- 
sis on phase-dependent mechanical force. Since we prove 
that this force emerges from charging effect, we will start 
with a detailed discussion of electrostatic energy in the 
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setup, and express the forces and superconducting cur- 
rent in terms of the induced charge. By doing so, we as- 
sume that the typical time of charge equilibration is much 
shorter than the typical time scale uJq^ of the mechanical 
motion. We compare the electrostatic phase-dependent 
force with Lorentz force proposed in 20 . We derive the 
coupled equations of motion governing the Josephson and 
mechanical dynamics and identify the dominant source of 
non-linear behavior. 



A. Electrostatic energy 

The junctions connecting the middle of the nanotubc 
to the leads have intermediate resistance so that the mid- 
dle of the nanotube forms a Coulomb island that is nei- 
ther isolated from nor ideally connected to the leads. 
While this situation is difficult to quantify from a mi- 
croscopic calculation, it can be completely analyzed at 
phenomenological level. The consideration below is just 
a case of elementary non-linear electrostatics and is sim- 
ilar to the discussion in |19j . However, it contains some 
important and less obvious details so we choose to present 
it at a comprehensive level. 

To start with, let us assume that the capacitance to 
the gate is vanishingly small while Vg is diverging such 
that the charge induced to the Coulomb island by the 
gate q — CgVg tends to a constant limit. A part of 
the ground-state energy of the setup, Ec{q) does de- 
pend on q. General Coulomb-blockade considerations 
[22] imply that this part is a (quasi) periodic function 
of q with a period of 2e. In the limit of full isola- 
tion, for instance, this energy is piecewise parabolic, 
Ec{q) = Ec'TXiini^ (N — q/2e)'^, N being an integer num- 
ber of extra Cooper pairs stored in the island. In general, 
it is a smooth function of q and may depend on the su- 
perconducting phase difference ip and, in principle, on 
mechanical displacement y. This energy results in a non- 
zero electrostatic potential difference between the island 
and the leads, V{q) = —dEc{q)/dq. 

Let us now turn to finite Cg and therefore finite Vg 
that is the potential difference between the leads and the 
gate electrode. Since this is not the potential difference 
between the island and the gate anymore, the induced 
charge q is not equal to C'gVg. Rather, it is determined 
from the voltage division in a nonlinear capacitance net- 
work, or, equivalently, from the minimization of the total 
electrostatic energy with respect to q, 



E = min, Ec{q) 



2C, 



-QVg 



(6) 



The charge is then found from the condition of the min- 



CgVg 



dEM 



(7) 



There are two implicit dependences in this equation that 
distinguish it from pure electrostatics, and that we make 



explicit now. First of all, the electrostatic energy depends 
on the mechanical displacement of the nanotube. Geo- 
metric considerations suggest that this dependence can 



be ascribed to C„ 



Indeed, the modification of capaci- 



tance to the gate is linear in y, Cg — )■ Cg -I- -j^y while 
the modification of E^ is expected to be oc y^. Second, 
the electrostatic energy depends on the superconducting 
phase difference: indeed, the Josephson energy is just the 
phase-dependent part oi E^, Ec{q,(p) = Ec{q) + Ej{q,(p). 
The electrostatic charge q depends both on displacement 
y and on superconducting phase ip. 

To single out these contributions, we assume that i. the 
voltage between the middle of the nanotube is smaller 
than the gate voltage, dEc/dq ^ Vg, this is fulfilled if 
the induced charge q^ e, i.e., in any practical setup; ii. 
the mechanical displacement is small in comparison with 

the distance to the gate, y <C Cg (^^^^ — Eg (in our 

experiments Lg c^i L ^ y). With this, we linearize Eq. 
([t]) with respect to the Josephson energy and mechanical 
displacement to arrive at (go = CgVg) 



(8) 



The first term is the common expression for the gate- 
induced charge while the second and the third are the 
corrections of interest. At the moment , we keep q in 
the argument of Ej, although q ~ qo- The point is that 
the Josephson energy is sensitive to variations of q of 
the order of e, and {q — qo) can in principle be of this 
order. Since y <^ Lg, we may disregard the possible y 
dependence of dCg/dy. 



B. Forces 

The mechanical resonator is affected by the electro- 
static force F = —dE/dy: 



F : 



dCg q^ 
dy 2Cl 



(9) 



It is convenient to distinguish three separate contribu- 
tions to the total force: the static force, the gate driving 
force and the phase-dependent Josephson force. 

The static force is produced by the d.c. gate voltage. 
Its magnitude is given by F^tatic = {dCg/dy)V^Q/2, cor- 
responding to the first dominating term in ([8|.The effect 
of the static force is to pull on the resonator, thereby 
tuning its resonance frequency [T3]. Since it is stationary, 
it does not excite the oscillations. 

The a.c. gate driving force arises due to the a.c. 
modulation of the gate voltage and is given by Eg = 

idCg/dy)Vg0Vg ~ {qO / c) (cVg / L g) . 

The phase-dependent Josephson force, not discussed in 
previous literature, comes about the product of the first 
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and third term in dSl) 



dCg dEj{q,Lp) 



dy 



^gO- 



dq 



(10) 



In fact it is similar to the gate driving force, with Vg 
replaced by the phase-dependent voltage arising in the 
capacitive network, dEj (g, tp) / dq. In contrast to the gate 
driving force, the time dependence of the Joscphson force 
is determined by the phase dynamics rather than the 
external modulation of the gate voltage. 

The scale of the Josephson force is Fj ~ {q/e){Ej / Lg), 
where we assume dEj/dq ~ Ej/e, which is true for inter- 
mediate contact conductances ~ Gq. It can be compared 
to the scale of a.c. gate driving force, Fj/Fg ~ Ej/eVg. 
For sufficiently low a.c. driving amplitude eVg <ti Ej the 
Josephson force dominates Eg <C Fj. 



C. Lorentz force 

The Josephson force explained above arises from the 
combined effect of charge sensitivity of the Josephson 
coupling and the capacitive coupling to the gate elec- 
trode. The alternative mechanism of generating a (p- 
dependent force was recently proposed by G. Sonne et 
al [50]. This force is of Lorentz type arising from the 
interaction of the (ys-dependent superconducting current 
with an external magnetic field B applied in the perpen- 
dicular direction (for the setup of Figure [l| along the 
2-axis.) This mechanism does not require the presence 
of a gate. 

Let us compare the Lorentz force and the electrostatic 
Josephson force. The Lorentz force is in y-direction, 
that is, perpendicular to both B and the supcrcurrcnt, 
Fb = L\B\Is. It is natural to express Fj in terms of 
the electric field \E\ = {dCg/dy)Vg/Cg produced by the 
gate electrode. For estimates, we assume L ~ Lg and 
{dCg/dy)L/Cg ~ L/Lg ~ 1. This yields 



Fb c\B\ 



F, 



\E\ 



(11) 



where c ~ 3 x 10® m/s is the speed of light and a = 
e^/4:TT€ohc ~ 1/137 is the fine structure constant. 

Typical magnetic fields used in experiments are |_B| ~ 1 
T. They are limited from above by the critical fields of 
the superconducting leads. The typical electric fields are 
\E\ ~ 10^ V/m. This corresponds to a potential drop of 
Vg ^ 10 V over a distance of Lg ~ 0.5 ^m. For these 
values Fb/Fj ~ 10a <C 1 suggesting that the Josephson 
force dominates. Therefore in the rest of the paper we 
will disregard the Lorentz force. 

If one imagines a ballistic nanotube, the Josephson 
coupling is not affected by the induced charge. In this 
case the Lorentz force would be the only superconduct- 
ing phase-dependent driving mechanism. However, the 
ideally ballistic nanotubes have not been realized exper- 
imentally. 



D. Analysis of non-linearities 

Let us bring together three coupled equations govern- 
ing the dynamics of the setup 



+ 1- 



Ze{uj) 2e Ze{u}) 



(^^(-7,^))^ , (12) 



y + ry + uj^y + ay' - f3y^' = M 



-idCg q' 
dy 2C2 



90 dCg dEj 



(13) 



(14) 



The first equation describes the dynamics of supercon- 
ducting phase difference ip{t) and is obtained by applying 
Kirchhoff's laws to the circuit. The second equation is 
for the mechanical displacement y{t) where we substitute 
the electrostatic force (|9|. The induced charge q enters 
both equations, and at the same time is defined by the 
third equation, that is, its value depends both on ip and 
y. Therefore the equations are coupled. 

We wish to simplify these equations under experimen- 
tally relevant assumptions. For this, we shall analyze 
the relative importance of different non-linearities in the 
coupling. There are natural non-linearity scales for all 
three variables, ip ~ 2tt, q ~ e, y ~ yc — -y/woF/ /?'. 
This could change if the coupling is sufficiently strong. 
For instance, the displacement may cause the variation 
of phase that is subject to Josephson non-linearity. The 
resulting variation of phase would produce the non-linear 
variation of q, this will result in non-linear feedback on 
y. This could in principle cause a non-linear scale of y to 
be smaller than j/c- So first of all we shall quantify the 
coupling between electrical and mechanical variables by 
comparing the non-linear terms in the mechanical force 
resulting from the coupling with those coming from the 
intrinsic non-linearities characterized by a and /3. 

The conclusion of this Subsection is that the mechani- 
cal non-linearity is the dominant one. We prove this with 
a rather involved reasoning given below. 

For the estimations, it is convenient to introduce the 
following dimensionless parameters: 



A. 



C. 



9 dq^ 

hujQ 



Ej_ 
Ec 



dq 



Ec 



- c. 



For estimations, we assume that Aj,Bj are either small 
or of the order of 1. This assumption is valid for Aj] 
it compares the Josephson energy to the charging en- 
ergy under conditions of well-developed Coulomb block- 
ade. The parameter Bj is a coefficient of Josephson feed- 
back at high frequency and depends on the details of the 
external circuit via the impedance Z^. Unless a special 
effort is made to increase the circuit impedance at high 
frequency, Bj will not be big. 

Given a variation of displacement dy we estimate the 
linear responses of the charge Sq and the superconducting 
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phase 5tf on 5y using Eqs. (14 1 and (12) as 



that are quadratic and cubic in 5y. 



6q{l 

6^ {I 



A,+A,B,) = B, 



Cg dy 
qo 1 dCg 
e Cg dy 



B 



<lo5y_ 
e L„ 



(15) 



(16) 



Assuming Aj ,Bj < 1 the hnear responses can be esti- 
mated as 



Sy 



qo5y_ 
e Ln 



We use this to find a scale of 5yc for which the responses 
of charge and superconducting phase may become com- 
parable with the scales of their intrinsic non-hnearities 
(5g ~ e and 5ip ~ 2tt. We conclude that 6yc — Lg{e/qo) 



Comparison with the scale yc 
non-hnearity yields 



Sy. 



r I \fQ of the mechanical 



2/c go r 



(17) 



Two last factors in this expression are big, while the first 
one can be small. We estimate the biggest go from the 
condition that cjq is changed significantly by applying the 
gate voltage, that is, the stationary displacement yo — 
This yields 



e 



(18) 



a being atomic scale, go — lO^e for our devices. The fact 
that the intrinsic scale of the non-linearity is smaller sig- 
nals that the intrinsic non-linearity of the resonator dom- 
inates the non-linearity arising due to coupling. With 
this, we estimate the first two factors as (e/go)(L/r) ~ 
{aL/r^Y. This is ~ 10 for our geometries and we con- 
clude that Sy/yc^ 1 for any Q > 1. 

This proves that the dynamics of charge and phase is 
linear in y provided our estimations of mechanical non- 
linearities a, f3 hold. We still need to show that the cou- 
pling to Josephson junction does not change these non- 
linearities significantly. 

So we estimate the quadratic and cubic non-linearities 
of the mechanical force due to coupling. First we find the 
quadratic and cubic variations of charge with respect to 
displacement using Eq. (14 1. 



go Sy 



Sl'''-Cg^{Sqrc.eA,(^ 



e Lr 



(19) 
(20) 



We can now estimate the terms in the mechanical force 



- dy 2Cl ' 




^stationary 




dy 2C2 \ go go go / 



t at i on ar y j 



(?) 



(21) 



(22) 



We compare these terms with the intrinsic non- 
linearities. The second order term SF'-'^'> needs to be 
compared with the mechanical quadratic non-linearity 
MaSy'^. Assuming the static displacement of the order 
of CNT radius, yo ~ r, we estimate Ma c:^ M/3/r ~ 



MaSy'^ 



< 1 



Here 



(23) 



1, 



we use the estimation {qo/e){r/Lg) <^ 
{qo/e){r / Lg) ~ 0.1 for typical CNT geometries, (see Eq. 
(18l). The third order term SF^^") needs be compared 



with the third-order non-linearity Mj35y . This yields 



MpSy^ 
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(24) 



To summarize, we proved that the non-linear scales 
correspond to 93 ~ 27r, q ~ e, y ~ yc — \J loqT / (3' and 
that for a CNT resonator the intrinsic mechanical non- 
linearities dominate the non-linearities arising from cou- 
pling. This permits a simplification of the dynamical 
equations. We may linearize the terms describing the 
coupling of mechanical displacement and electricity, thus 
separating Josephson and mechanical non-linearities. 



E. Workflow 

This sets the following workflow. 

• At given a.c. and d.c. bias and gate voltages we 
solve for Josephson dynamics neglecting the me- 
chanical coupling and setting q = qo{t). We flnd 
I{t) and ip{t). Using these, we compute the Joseph- 
son force Fj given by Eq. (10). 



We solve the non-linear mechanical equation 



M{y + Ty + uj^y + ay^ - /Jy^) 



Fst + Fg + F, 



(25) 



to find y{t). We are mostly interested in a part 
that oscillates with frequency ~ loq. This may be 
excited by both Fj and Fg. 



8 



• In most cases, we are interested in a (d.c.) current 
response on the mechanical motion, the mechani- 
cal response. It arises due to direct modulation of 
charge by the mechanical displacement, 



(26) 



in the first order in y{t). We will mostly concentrate 
on the situation when the displacement oscillates 
at the resonant frequency while a d.c. component 
of /mil is of interest. The d.c. signal them arises 
from the rectification of y{t) by an oscillating part 
of dl{(t>)/d{q/e), that we call the detecting current. 

If we can neglect the feedback in Josephson dy- 
namics, we are done, since the response is given di- 
rectly by /mh- Otherwise, we linearize the Joseph- 
son dynamics to determine the response of the su- 
perconducting phase on the mechanical displace- 
ment found, (^mh(i) 



2e 



dt'dt"Z{t',t"%,Y,{t") 



Here the kernel Z(i, t') represents the combined lin- 
ear impedance of the external circuit and the junc- 
tion. 

• Taking this into account, we obtain the current 
response sought: 



_ ~ 2e d'^Ej 

/mh — /mh + -^-g^fmh 



(27) 



The first term is the direct modulation of the cur- 
rent by the charge induced by the mechanical dis- 
placement while the second one is a feedback of 
Josephson junction by means of (pmh- 

The mechanical response is thus typically a small cor- 
rection to the maximum superconducting current. We 
can estimate it at maximum taking the 



/m 

T 



qoVc 

e L 



QeL 



(We remind that {qo/e){r/L) ~ 0.1, Q ~ 10^ for our 
devices). Perhaps unexpectedly, the typical response be- 
comes smaller upon increasing the quality factor. The 
reason for this is clear: the maximum displacement be- 
comes smaller. However, the large Q results in sharp 
frequency dependence of the response making it easier to 
identify. We thus concentrate on this dependence. 



F. Parameters 

Let us specify the values of parameters employed. We 
choose these values such that they closely match those of 



a typical experimental realization of our setup. 

Junction critical current, Ic : 1.0 x 10^^ A; (28) 

Josephson energy, Ej : 2.1 x 10"*^ eV; 

Static gate voltage, Vgo : 1.0 V; 

Static charge on the resonator , qo/e : 100; 

Josephson force, Fj : 1.1 x 10^^^ N; 

Resonator length and distance to gate, L = Lg : 0.3 /zm; 

Resonator mass, M : 4.1 x 10^^^ kg; 

Resonance frequency, wo/27r : 0.84 GHz; 

Quality factor, Q : 1.4 x 10^; 

Quadratic non-linearity, a : 5.5 GHz^nm^^; 

Cubic non-linearity, /3 : 5.5 GHz^nm^^; 

Scale of maximum displacement, yc : 3.2 pm; 



Mechanical force scale, Fr : 1.2 x 10 
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IV. PHASE BIAS 

Let us start our considerations with the junction bi- 
ased with a time-independent phase (p: such bias condi- 
tion can be achieved by embedding the junction into a 
superconducting loop. Unfortunately, our present exper- 
imental setup does not allow measurements under these 
bias conditions. We present the theoretical results in 
hope that they will be useful for future experiments. 

The simplest experimental signature of Josephson force 
under phase bias conditions is the phase-dependent shift 
of the resonant frequency. The mechanism of this shift in 
our situation is the mechanical non-linearity: the static 
Josephson force tightens or looses the nanotube resulting 
in the frequency change. The frequency response on the 
static force in our model reads 



dujQ 
~dF 



Wo 



MujI 



^static 

SO the phase-dependent frequency shift reads 



(29) 



a dCg dEj{q,ip) 



Muj^ dy 



dq 



E, 



OJQ 



Wo 



static 



eVr 



COs((/3) 



(30) 



The shift is clearly observable provided it exceeds the 
broadening T. The estimation gives 



— — ^ Q— ^cos^ . 
r eVgo 



(31) 



For the parameter set ( 28 ) , maximum value of the shift 



= 2.8, 



V r / 

this is, the shift is clearly observable. 



(32) 
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Let us consider an example of a mechanically-induced 
response under conditions of the phase bias. To excite 
mechanical oscillations, we apply an additional a.c. volt- 
age to the gate that oscillates at the frequency 51 close to 
the resonant frequency ujq- 



(33) 



Assuming Vg to be sufhciently small to provide a linear 
response of the displacement, we obtain the following ex- 
pression for the resonant part of the displacement: 



y 

y(.t) 



-1 



2Mwo i'i(p) + iT/2 
-lye +y e ) 



dC, 
dy 



ffOfff , (34) 
(35) 



Here, v{ip) = 57 — wq ~ Awq ((/?) <C wq is the detuning that 
includes the phase dependent shift of the resonance fre- 
quency discussed above. Owing to the mechanical non- 
linearity, the oscillating displacement produces a station- 
ary displacement y — a\y\'^ /luq. This induces a station- 
ary charge that affects the d.c. superconducting current 
at constant phase bias. Rather remarkably, this effect is 
related to the phase-dependent shift discussed above. In- 
deed, both are proportional to charge-dependent part of 
the Josephson energy and to the non-linearity coefficient 
a. The resulting current response reads 



-2e— {^ojoiif)) lyp(AfwoM), 



(36) 



the contribution to the current being of the order of 
eAwo if we measure quantum fluctuations of displace- 
ment ^JKJMjjJq. The dependence on frequency of the 
a.c. modulation is a Lorentzian one, as it is frequently 
expected (Fig. |3]), the Lorentzian center being shifted 
with changing the phase. 
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FIG. 3: (Color online.) The phase-dependent frequency shift 
for the case of weak driving. The curves give the hnear re- 
sponse of d.c. current in units la ~ 2e(Aa;o)max|j'P(Afa;o/ft) 
as a function of frequency detuning for a set of phase bias 
values: from the lowermost to the uppermost curve the phase 
changes from = J to = ^'-l^, with interval 7r/8. The 
curves are offset for clarity. Dashed lines give the positions of 
zero. 



the linear regime. 



2 z/(^)2 -)-r2/4' 



(38) 



so that the signal changes sign at the resonance point. 
Fig. |4] illustrates the Fano-type dependence in the non- 
linear regime. Comparing expressions ( 38 1 and ( 36 1 we 
conclude that the Fano-shaped /„i dominates provided 
K)/^90 ^ ('Zo/e)"'^^^ — 10^'^, this is, deep in non-linear 
regime. 



D.C. VOLTAGE BIAS 



A. Fano-type response 

The above mechanism of response exploits the domi- 
nating mechanical non-linearity. It is proportional to y^ . 
Upon increase of the a.c. amplitude Vg the oscillating 
displacement saturates owing to the non-linearities. In 
this case, the dominating d.c. current signal can arise 
from the electrical non-linearity as a result of a mixing of 
the oscillating displacement y and the oscillating charge 
oc Vg. The resulting current is thus proportional to Vgy 
and may exceed the contribution oc y provided the latter 
saturates. 

The expression for this contribution reads 



2ed^Ej{q,ip)^ f dC, 



'-dip 



Ca 



dy 



Interestingly, it exemplifies a Fano-type dependence on 
the detuning that is quite different from a Lorentzian. In 



Let us turn to d.c. voltage bias. In this case, the 
superconducting phase is in first approximation a linear 
function of time, (p = ujjt, LUj = 2eV/h being the Joseph- 
son frequency which corresponds to the voltage V across 
the junction. In the same approximation, the current is a 
purely oscillatory function of time. The time-dependent 
current can be expanded into harmonics of the Josephson 
frequency. 



, , 2e dEj , , ^— \ f . , X 

l[t) = — ^— (q, = Wj<) = y J„ Sin(7T,u}jt). 



h d(p 



(39) 



A d.c. current emerges from the feedback on Josephson 
dynamics: oscillatory currents produce oscillatory cor- 
rections to the phase proportional to the impedances at 



frequencies 



(37) approximation in Z(uj), we arrive at 



Taking this into account in the first 



(40) 



10 




FIG. 4: (Color online.) An example of Fano-type frequency- 
dependence of the mechanical response (371. The curves cor- 
respond to the driving force values F/Fc = 0.2,0.6, 1 in the 
upper panel and F /Fc = 2, 6, 10 in the lower panel. The cur- 



rent is in units of 



loVgUc that amounts to ~ 10 ^/c 



dq^ dy 

for the parameter set chosen 



The above relation holds provided TLeZ ^ V/Ic- Non- 
perturbative treatment of Josephson dynamics is re- 
quired otherwise. 

Let us consider the mechanical effects. It is important 
to note that under the d.c. voltage bias the Josephson 
force also oscillates in time, 



Flit) 



dC dE,{q,v) 



dy 



dq 

dq 



cos{nujjt) 



(41) 



Ej^n being the harmonics of Josephson energy. Therefore, 
the force can efficiently excite the mechanical resonator 
provided nujj ~ ujq. Let us first concentrate on the case 
where the resonance frequency is matched by the first 
harmonics, LUj ~ ujq. The detuning is defined as ^ = 

UJj - UJQ. 

To start with, let us assume that the Josephson force 
is sufficiently weak so that the mechanical response is 



linear and thus given by Eq. ( 34 ) . The direct mechanical 



contribution to the d.c. Josephson current is obtained by 
averaging Eq. |26[ and reads 



dh fdC„ v. 



d{q/e) \ dy e 



Imy 



This can be cast to the form similar to (|40|) 



dh 



d{q/e) 



V 



(42) 



(43) 



where the current is replaced with detecting current 
dl/d{q/e) and the "mechanical impedance" Zmh{'-^) is 
defined as 



—IV + 1/2 

^ n^f dCg^VgV h 

e-^\dyej 2Mwo ' 



(44) 
(45) 



(Here, i/ = ujj — ujq). 

This form of the presentation of the mechanical re- 
sponse makes evident an analogy with Fiske steps [23] 
that are observed at voltages corresponding to resonant 
frequencies of an electrical impedance. This may be 
either an impedance of external circuit or an effective 
impedance that is essentially contributed Josephson in- 
ductance. 

To comprehend the scale of the response, we note 
first that for a sufficiently well-developed Coulomb block- 
ade /]_~ dli/d{q/e). Therefore, to compare the cur- 



rent (40 1 and the mechanical response, we need to com- 



pare and a typical environmental impedance. The 
latter can be estimated as the impedance of free space 
Zf ~ lO^Ohm. The typical mechanical impedance far 

from the resonance, Zm\ should be much smaller than 



that. Indeed, if we substitute the parameter set (|28|) into 

.(0) 



Eq. 



45 



we end up with Zm' = 0.7- 10~^ <C Zf. However, 
the "mechanical" impedance is enhanced by a factor of Q 
at the resonant frequency. With this, Z^^ > Zf and the 
current peak produced by the mechanical response should 
exceed the background current given and be clearly ob- 
servable. 

The voltage dependence of d.c. current response in 
linear regime is determined by KeZm and thus takes a 
Lorentzian shape with the half- width SV = V/Q. This 
assumes a noiseless voltage source. It is known |21j that 
the voltage noise suppresses the coherence of Josephson 
generation. For white noise spectrum of intensity Sv, the 
resulting line-width reads SVn — {2e/h)^Sv- Comparing 
the two, we conclude that the mechanical response will 
be broadened by the noise and essentially suppressed pro- 
vided SVn > SV, this is, Sv > {2e/H)ujj/Q. This gives a 
condition on the possibility of detection of the mechani- 
cal response that may be challenging to meet in practical 
circumstances. 

As mentioned, the Josephson force can be big enough 
to exceed Fc, this makes it relevant to address the non- 
linear response as well. We illustrate the non-linear re- 
sponse in Fig. [5) To produce the Figure, we took the 
charge-dependent part of the Josephson energy to be of 
the form E^^\q,ip) = i?cos((^) cos(7rq/e). Tuning q with 
the d.c. gate voltage tunes the Josephson force from to 
some maximum value. We chose E such that the maxi- 
mum force equals 10 Fc {E = 0.23 /iteV for the parameter 
set in use) and compute the response using Eq. ( 42 1 and 
Eq. ([5]) at a set of the values of q, or, equivalently, Fj. 
The response is Lorentzian at small forces, increases and 
develops a jump characteristic for bistability. It is in- 
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FIG. 5: (Color online.) Top left panel: The charge-dependent 
part of Josephson energy as function of the gate-induced 
charge q = CgVg. The crosses indicate the values of q that 
correspond to the values of Josephson force used in other pan- 
els. Top right panel: the mechanical response as function of 
detuning for relatively low values of the force F = Fj / Fc > 1, 
at which the response increases with increasing the force. 
Bottom panels: Frequency dependence at the force values 
F/Fc = 3, 5, 7, 10 (from the left to the right panel) where the 
response decreases with increasing the force. 



teresting to note that the response slowly decreases upon 
increasing Fj at Fj > 2Fc. This is because the response is 
proportional to Imy that quickly decreases at big driving 
forces. In this limit, Im oc F; 



B. Parametric excitation 

For the sake of completeness, let us mention the possi- 
bility of the resonant mechanical response at higher volt- 
ages by means of parametric excitation [T^]. Generally, 
parametric resonance in a non-linear oscillator is achieved 
by applying an a.c. driving force with frequency about 
a double of the resonant frequency, fl ~ 2u!o [12] . In our 
case, this is achieved by applying a d.c. bias voltage with 
ujj ~ 2ujo, so that the Josephson force oscillates at 2ujo 
and integer multiples of this frequency and thus provides 
the parametric driving required. 

The response of at resonant frequency emerges pro- 
vided the parametric driving force exceeds a certain 
threshold value, and, as in case of direct resonance, 
achieves values ~ yc- The point is that this threshold 
driving force is parametrically bigger than F^ Ft ~ \/Q- 
For our devices aX Q = 10^, the parametric excitation re- 
quires Ej that by a factor of 30 exceed the value from the 
parameter set and are not practical. This is why we do 
not explore the regime of parametric excitation in detail. 

Besides, the manifestation of the oscillating amplitude 
is not as straightforward as in the case of direct reso- 
nance. The contribution of displacement at Wq to the 
mechanical current response ( 26 1 oscillates at the same 
frequency and is not readily rectified to a d.c. current. 
Under our assumptions, the d.c. mechanical response 
is dominated by the displacement oscillating at 2uq and 
is by a factor of ^/Q smaller than the typical responses 
studied in this paper. 



VI. SHAPIRO STEPS AT RESONANT DRIVING 



A. Excitation by higher harmonics 

If we take higher harmonics of current-phase character- 
istic into account, we note that Josephson force emerges 
at a set of frequencies that are integer multipliers of ujj 
(Eq. (41)). This implies that the resonant mechanical 
response can be also observed in the vicinities of a set of 
discrete voltage values satisfying ujj = oJo/n^ this is, at 
lower voltages than the resonance described above. The 
response is computed along the same lines with replacing 
Ej i by Ej n- In linear regime, the response reads 



din 



d{q/e 



nReZmh {nujj 
V 



(46) 



(cf. Eq. 43 the factor n in the present expression is 
canceled by lower voltage V = {h/2e)u!Q/n) The response 
scales with the relative values of the harmonics and is in 
principle of the same order of magnitude for several low 
harmonics. Its dependence on voltage in the vicinity of 
the resonance is similar to that discussed above and does 
not have to be illustrated separately. 



From now on, we turn to the situation where the setup 
is a.c. driven at frequency fl. As discussed in Section 



II A in our setup this gives rise to two a.c. signals Vg (t) — 
VgCos{nt + x), Vb{t) = Vbcos(m). The effect of Vt is a 
formation of Shapiro steps [24] . 

A common approach to Shapiro steps takes into ac- 
count only the first harmonics of the current-phase re- 
lation and starts with the assumption that the time- 
dependent superconducting phase difference can be pre- 
sented as a sum of three terms 



ip{t) = (fii sin(rit) -f- LOjt + ifQ. 



(47) 



Here, the first term describes the a.c. driving {(pi > 0, 
fi — \Vb\/{2e/h)fl, the second term corresponds to a d.c. 
voltage V — LUj / (2e/h), and the third term is a lock-in 
phase important for further consideration. With this, 
sm{ip) can be presented as a sum over harmonics 



sm{ip) 



E 



Jm{'Pi)sm{nmt + (fio) (48) 



with n„i = mQ + ujj. Here, denotes the m-th Bessel 
function of the first kind. 
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Shapiro steps are formed at discrete values of d.c. volt- 
age \ujj\ = mfi. In this case, the time-dependent current 
I{t) = IcSm{ip{t)) has a d.c. component 



Making use of Eqs. (481 and (26 1, we express the me- 



/dc = -IcSgIl{V)Jm{ipi) Sinifio 



(49) 



Simplest assumption is ideal current bias at zero fre- 
quency and ideal voltage bias at frequencies ~ 51. In 
this case, the I-V curve of a.c. driven junction consists 
of a series of separate pieces. At each piece (Shapiro 
step) the voltage is locked to one of the discrete values. 
The current within each piece may vary from minimum 
values /_ to the maximum value /+ provided the bias 
current fits this interval. In this case, the actual value 
of the lock-in phase (pQ is set by the bias current. The 
extremal values I± = ±Ic\Jm{'Pi)\ are achieved at the 
lock-in phases given by 



(50) 



We in main follow this approach while admitting ex- 
treme simplifications it brings. The higher harmonics of 
current-phase relation and/or non-ideal voltage bias not 
only modify the relation between the current and lock-in 
phase: They also provide phase-locking at fractional ra- 
tios of ujj/fl |26j and formally at all rational values of this 
ratio. These fractional Shapiro steps are however more 
sensitive to noise than the integer ones and more likely 
to vanish. The I-V curves of our devices do show well- 
developed steps at integer values of ujj and only traces 
of phase-locking at intermediate values. For this reason, 
we do not consider fractional Shapiro steps in this paper 
and concentrate on integer ones where \ujj\ = mfl. 

It is advantageous to look at the mechanical response 
at Shapiro steps rather than at d.c. bias conditions. The 
external a.c. driving synchronizes Josephson oscillations. 
The inductive response present at Shapiro steps also re- 
duces significantly the voltage noise at low frequencies so 
that it does not broaden the resonant lines. 

In this Section, we will consider the mechanical re- 
sponse in the simplest situation of resonant driving where 
the driving frequency matches the resonant frequency, 

r2 ~ cjo- 



A. First step 

Let us first concentrate on the first Shapiro step, the 
one at voltage 2eV/h = ±51, that is the biggest in the 
limit of small driving voltages (pi ^ I, and determine the 
d.c. part of the response at the oscillating displacement 
y. To represent the results, we normalize y to the non- 
linearity scale yc and introduce a convenient current scale 



For the values of our parameter set, 

I ~h{qa/e){yjLg) = 1.1 x 10-=* 



chanical response in terms of the amplitudes y and Ij at 
the resonant frequency. 



-Re{j*y} 



(52) 



J2{pi)e'V'> - Jo(^i)e-^'^° 
Jo((^i)e*'^« - J2{vi)e-'^'> 



if > 
if F < 



(53) 



This displacement is a response on the force at resonant 
frequency which is a sum of Josephson force and gate 
force. The timc-dcpcndent Josephson force is expanded 
in harmonics in the form 

F,(^) = F,cos(^(t)), (54) 

OO 



m— — OO 



dy 



go- 



dq 



e La 



Its amplitude at resonant frequency is contributed by the 
terms m — Q,2 and reads 



Fj = Ff; 



J2(¥'i)e"^'' + Jo(<^i)e-^'^") ifF>0 
Jo(¥>i)e''^'' + J2(v?i)e-*'^") ifF<0 



(55) 
(56) 



Let us discuss first the relative scale of the gate force in 
comparison with the Josephson force. It may seem there 
is none, and varying the a.c. gate voltage Vg one can 
achieve any ratio ~ Vg/Ej between the forces. However, 
we should take into account the fact that in our setups 
a.c. driving also induces an appreciable bias voltage Vf,. 
If the oscillating phase (pi produced by this voltage be- 
comes large Shapiro steps can hardly be observed. It is 
in general reasonable to expect Vg ~ H- In this case, 
(y9i ~ 1 corresponds to Fg/Fj ~ huo/Ej. The latter ratio 
is typically 10^^ in our setups (for our parameter set it is 
huJo/Ej — 2.7 X 10^^). This implies that typically we can 
disregard a.c. gate force in comparison to the Josephson 
force. We will analyze this case first and consider the 
effect of the gate force in the end of the Subsection. 

With this, the mechanical response is given by 



IF 

/mh = —RciffR} = 

i'c 

sgnV^(j2(^i)- Jl(vPi)) Im(i?) 

- Jo((^i)J2(¥'i)sin(2v7o)Re{i?}) (57) 



IF 



(51) Here, R = R{v / {T /2),{Fj / Fc)\f\) defined by Eq. ^ 
gives the non-linear mechanical response. The expression 
is naturally separated onto two terms. The first term is 
proportional to Im(i?) and therefore exhibits a Lorentz- 
like dependence on frequency. It does not depend on the 
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lock-in phase and can be regarded as a shift in the cur- 
rent. Owing to the shift, the maximum and minimum 
currents at the step are no more opposite: the me- 
chanical effect breaks the symmetry of the Shapiro step. 
The shift is however opposite for opposite voltages. The 
origin of the shift may be traced to the Fiske response 
(Eq. 43 ) formed at ujj ~ ujq in the absence of the a.c. 
driving. Indeed, in the limit of vanishing ipi the first 
term in the mechanical response does not vanish: rather, 
it approaches the expression (43). So it looks like the 



Fiske response persists also for well-developed Shapiro 
steps and contributes to the current at the step. This 
suggest perhaps the easiest way to observe and identify 
the mechanical response: measure maximal and minimal 
currents at a step as function of the a.c. frequency. In 
the rest of the paper we thus mainly concentrate on the 
modification of extremum currents. 

The second term in Eq. (57) cannot however be ob- 
served in this way. In ideal current bias conditions at 
low frequency, the second term in the current response 
fact amounts to a shift of the lock-in phase. Indeed, 
since the current at the step as function of ipo reads as 
I {if) — — sgn(F) Ji(<y9i), the second term can be seen as a 
modification of the lock-in phase at constant bias current 
which does not depend on this current. 



(A(^, 



0)mh 



-sgn(y) 



IF Jo(yi)J2(yi 



.Re(i?) (58) 



This response is of Fano-type. Since such shift of the 
phase does not modify the values of the current ex- 
trema, the effect cannot be observed in the course of 
two-terminal electrical measurement in our setup. The 
shift of the lock-in phase can be however revealed if the 
Josephson junction under consideration is a part of a 
SQUID, or with the aid of lock-in measurement at non- 
resonant a.c. frequency. 

With respect to this, we ought to mention yet another 
effect of Josephson force manifesting itself in the me- 
chanical response considered. In fact, the situation at a 
Shapiro step is similar to the phase bias conditions con- 
sidered in Section [TV] with lock- in phase playing the role 
of ip. We thus expect (/Jo-dependent shift of the resonance 
frequency. The static Josephson force at a Shapiro step 
is given by 

=i^sgn(l/)Ji(¥.)cos(¥>o) (59) 

The frequency shift caused by this force thus reads 

Aa;o((po) = sgn(F)(Awo)max-^i(¥') cos(<^o)- (60) 

Here (Acjo)inax is a maximum frequency shift in the ab- 
sence of the a.c. driving, given by Eq. (30). The fre- 
quency shift vanishes at extremum points (/?g and there- 
fore cannot be observed by measuring the extrema of the 
current. 

We illustrate the mechanical response in Fig. [6j In 
this Figure as well in all subsequent Figures except Fig. 




FIG. 6: (Color online.) The mechanical response at resonant 
driving ~ ojo and at the first Shapiro step Vb = 2en/ft 
versus the oscillating phase fi. The first (upper) plot gives 
the maximum current at the step. The second and fourth 
plots give the mechanical response defined as the extremum of 
the modification of this maximum current over the frequencies 
near the resonance, at maximum Josephson forces F = Fc and 
F = 50Fc, respectively. The actual amplitudes of the resonant 
Josephson forces are given at the lower plots, respectively 
third and fifth. 



[Tj we concentrate on the modification of the maximum 
current on the step. Instead of presenting the (rather 
trivial Lorentz-like) frequency dependence of this modi- 
fication, we give the extremum of this modification over 
the frequency range and plot it versus fi . The extremum 
is proportional to the maximum of Im(i?) over the fre- 
quency. We shall note that the dependence of this max- 
imum on the force is rather specific one: it is a constant 
until the bistability threshold F = 1.24i^c, has a cusp at 
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this value of force, and decreases monotonously at higher 
forces. This accounts for rather strange appearance of 
the response curves. If the Josephson force is smaller 
than the bistability threshold, they coincide with the lin- 
ear response given by the dotted curves. Otherwise, the 
response is smaller than linear one and exhibits kinks. 

The so-defined maximum response is plotted in Fig. 
[g] for two values of F, those correspond to slightly and 
strongly non-linear regime, respectively. In both cases, 
the response vanishes when the width of Shapiro step 
reaches maximum, or becomes zero (except (fi = 0). 
In slightly non-linear regime, the response reaches maxi- 
mum value at ipo. Upon increasing tpi, it exhibits Bessel- 
like oscillations. The envelop of these oscillations shrinks 
with increasing ipi. This shrinking is much faster than 
that for either step width or Josephson force. In strongly 
non-linear regime, the amplitude of the response is deter- 
mined by competition of two factors: it is increased by 
the bigger value of the detecting current, and decreased 
owing to smaller imaginary part of oscillating displace- 
ment at higher Josephson forces. 

Let us turn to the effect of the gate force. The Joseph- 
son force of the kind considered can change sign and 
therefore be tuned to zero by tuning q. In the vicinity of 
this particular q, the gate force should compete with the 
Josephson one and eventually dominate. The full ampli- 
tude of the force at the resonant frequency then reads 



F^Ff + FgeM-^x), 



(61) 



the frequency shift x between the bias and the gate volt- 
age being a relevant parameter. 

The mechanical response is given by Eq. ( 57 1 where R 



depends on the full force plus an addition proportional 

to Fg, 

4'h = ^M/exp(-zx)i?} = 

IF 

— (-(Jo(<^i) sin((/7o + x) + Mvi) sin{ifo - x)) Re(i?) 

(Jo((/9i) coa{ipo + x)- 'hifi) cos(.^o - x))Im(-R)) 

(62) 

The last equation holds for V > Q. The expression for 
1^ < is obtained by interchanging Jq a-nd J2. Evaluat- 
ing this at the extremum points of lock-in phase, (f^ , we 
obtain 



(9)± 



Tsgn(FJi ((/?)) 



IF„ 



(Joifi) + >/2(¥'i))Ile{i?exp(-ix)} 



(63) 



Therefore, the contribution of the gate force to extremum 
currents is not like a shifts: rather, it modifies the width 
of the step. These terms are even in voltage and display 
a mixture of Fano-type and Lorentzian-type response as 
function of frequency, this being tuned by the phase x- 

To illustrate a rather complex interplay of Josephson 
and we plot in Fig. Q the frequency dependence of the 




FIG. 7: (Color online.) The effect of gate force. For all 
plots, the gate force is fixed to Fg = F^. In the plots from 
top to bottom the maximum Josephson force F assumes the 
values F/Fc — —5, —1, 0, 1, 5. We choose x = ^^nd ipi — 1. 
Dashed fines in the plots for I~ give values opposite to the 
corresponding 7^, to stress the symmetry or asymmetry of 
the response. 
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mechanical response for a constant Vg and a set of values 
of Fj that pass zero. The plots show the modifications 
of extremum currents These modifications are the 
same for the Josephson force contribution and opposite 
for the gate force contribution. Besides, the frequency de- 
pendence is Fano-like for the gate force contribution and 
Lorenz-like for the Josephson force contribution. In the 
central plot, the Josephson force contribution is absent, 
the modifications of are opposite, and the frequency 
dependence is Fano-like. Upon increasing the Josephson 
force, these features are transformed into the opposite 
ones. The plots are symmetric upon simultaneous change 
of signs of the current and the Josephson force. 



B. Higher steps 

At the same conditions of the resonant driving, we ana- 
lyze the mechanical response at other Shapiro steps m > 
1, those correspond to higher voltages \V\ = ■mhujo/2e. 
Both the amplitudes of the detecting current and the 
Josephson force display a complex dependence on the 
step number m and the oscillating phase (fi. They are 
given by 



j = -i ( J-l+me ^° - Ji+rrie 



1+rn 



Jl 



-tipo 



where the dependence on the sign of the voltage is in- 
corporated into ffi = — sgn(V^)TO. We consider only the 
situation when Josephson force dominates. With this, we 
obtain a relation similar to Eq. (57): 



mil 



IF 

~f: 



ReiffR) 



IF 



{sgaV {Jl^,{ipi) - Jl+,{vi)) Im(i?) 
- Jm_i(<^i)Jm+i((pi)sin(2(^o)R.e(i?)) (64) 

As in the first step, the response consists of two terms. 
The first one gives a shift in the current, and gives a 
modification of the maximum and minimum currents at 
the step, this is to be measured. As in the previous case, 
the shift is odd in voltage. However, its (^i-dependence 
is quite rather distinct. 

The measuring of the mechanical response at higher 
steps is important to check the consistency of results and 
thereby unambiguously identify the mechanism of the re- 
sponse. The characteristic dependences on ipi make the 
identification easy. 

We illustrate the response for higher steps in Fig. [8] 
(second step) and Fig. VI B (fifth step). In both cases. 



the response correlates with the Shapiro step width given 
in the upper plots: it vanishes when the width achieves a 
maximum or becomes zero. In distinction from the first 
step, the responses vanish at vanishing ipi. Their typical 
values are of the same order. However, the envelopes of 
the responses decrease rather slow with increasing ipi. 




FIG. 8: (Color online.) The maximum current at the Shapiro 
step, the mechanical response and the force versus <fi for F = 
Fc and F = 50Fc at the second Shapiro step Vb = 4en/fi. 



VII. SHAPIRO STEPS AT NON-RESONANT 
DRIVING 



In the previous Section, we concentrate on the case 
when the driving frequency matches the resonant fre- 
quency of the mechanical oscillator. It is not a neces- 
sary condition for an efficient excitation of the resonant 
mode. The Josephson dynamics at Shapiro steps are es- 
sentially non-linear. As a consequence, the spectrum of 
current oscillations contain all higher harmonics nfl of 
the driving frequency fl. The same pertains the Joseph- 



16 




FIG. 9: (Color online.) The same as in Fig. [s] at the fifth 
Shapiro step Vb = lQeQ./h. 



son force. Therefore, the resonator can be efficiently ex- 
cited for 51 = loq/N, N > 1 being an integer number. 
At any given N, the resonant conditions are achieved 
for any Shapiro step number m, and thus for vohages 
2eV/h = ujj ^ {m/N)LOQ. 

These non-resonant driving conditions are advanta- 
geous for observation of the Josephson force since the 
a.c. gate voltage force is not in the resonance, does not 
cause any appreciable displacement and therefore does 
not mask the effect of the Josephson force. In this short 
Section, we will thus concentrate on the case of the non- 
resonant driving fl = ujq/N. 



The amplitudes of the detecting current and the 
Josephson force depended not only on the step number 
m and the oscillating phase (pi , but also on N. They are 
given by 

where the dependence on the sign of the voltage is again 
incorporated into m = — sgn(y)m. 

Since the gate force is absent, the response is given by 



a relation similar to Eq. 57 that contains the Josephson 
force only: 



IF 



ReiffR) 



IF 



{agnV (J^_Ar((pi) - Jl,+Nifi)) Im(^) 

- ^i-Ar('(3l)Jm+Af(¥'l)sin(2(po)Re(i?)) . 



(65) 



It is again separated onto two terms discussed above, 
only the first term being responsible for the modification 
of the extremum currents of the Shapiro steps. 

We illustrate the dependences on ipi in Fig. 10 (for 
N = 2) and Fig. [Tl] (for TV = 3) for the first three steps 
with m = 1,2,3. The vertical arrangement of the plots 
is the same as in the previous Figures except we chose a 
single value of the maximum Josephson force F = 50Fc 
that brings us deep into the non-linear regime. In this 
regime, the response is of the same order of magnitude for 
all steps and ratios iV, while retaining unique m, N spe- 
cific dependence on ip that can be used for identification 
of the effect and the characterization of the Josephson 
force. 



VIII. CONCLUSIONS 

In conclusion, we have studied Josephson junction dy- 
namics affected by excitation of a mechanical resonator. 
We have demonstrated that the mechanical oscillations 
can be rectified giving rise to an additional d.c. current 
that can be used for detection. The mechanical response 
is proportional to the oscillation amplitude, and is es- 
timated as /mh ^ Ic{qo/e){y/L) ~ IQ-^Idy/yc)- The 
resonator can be driven by the a.c. voltage applied to 
the gate electrode as well as an additional mechanical 
force, termed the Josephson force, that depends on the 
superconducting phase difference at the junction. We es- 
timate the Josephson force as Fj ~ {qo/e){Ej / Lg) and 
show that it is sufficiently strong to drive the mechani- 
cal resonator into the non-linear regime. We also show 
that it is typically larger than magneto- induced force pro- 
posed in pO] . 

We have presented a general and detailed analysis of 
the coupling between electrical and mechanical degrees of 
freedom, discussing the competing non-linearities. This 
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FIG. 10: (Color online.) The mechanical response at the non-resonant driving. Here, the a.c. driving frequency is f2 ~ tJo/2, 
corresponding to A'^ = 2. From left to right the three columns correspond to Shapiro steps m = 1, 2, 3. Plotted are the maximum 
current at the Shapiro step, the mechanical response and the amplitude of the force. The maximum of the Josephson force was 
set to F = 50Fc for all plots. 



^ 



fa 




FIG. 11: (Color online.) The mechanical response at the non-resonant driving for A'^ = 3. Except this, all other parameters 
are the same as in Fig. [To] 



analysis is applied to a Josephson device with a sus- 
pended CNT resonator, where we show that the intrinsic 
non-linearity scales dominate those arising from the cou- 
pling. 

We have provided analytical formulas for the response 
of the device to mechanical excitations in a wide inter- 
val of the excitation strengths and for various biasing 
schemes. We discuss distinct frequency dependencies, 
Lorentz and Fano-like, of the mechanical response both 
for linear and non-linear regimes and show how these 
arise based on the nature of the resonant mechanical 
force. In the case of a phase biased junction we show 



that the resonant frequency of the mechanical mode ac- 
quires a measurable phase-dependent shift (see Fig. [s]). 

We have discussed conditions of detecting the en- 
hanced mechanical response arising when the Josephson 
frequency matches the resonance frequency of the me- 
chanical mode. We reasoned that the regime of Shapiro 
steps is advantageous, since the fluctuations of the volt- 
age drop over the junction are suppressed. We provided 
expressions for the mechanical response in the regime of 
Shapiro steps and demonstrated that it manifests as mod- 
ifications of the extrema of the steps. We show that the 
mechanical mode can be efficiently excited not only by 
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resonant a.c. signals, but also by a.c. signals with fre- 
quencies close to an integer fraction of the mechanical res- 
onance frequency. Our preliminary experimental results 
confirm this behavior; these will be reported elsewhere. 
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